Explosive synchronization transitions in complex neural networks 
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It has been recently reported that explosive synchronization transitions can take place in net- 
works of phase oscillators [Gomez- Gardefies et al. Phys. Rev. Letts. 106, 128701 (2011)] and chaotic 
oscillators [Leyva et al. Phys. Rev. Letts. 108, 168702 (2012)]. Here, we investigate the effect of a 
microscopic correlation between the dynamics and the interacting topology of coupled FitzHugh- 
Nagumo oscillators on phase synchronization transition in Barabasi- Albert (BA) scale-free networks 
and Erdos-Renyi (ER) random networks. We show that, if the width of distribution of natural fre- 
quencies of the oscillations is larger than a threshold value, a strong hysteresis loop arises in the 
synchronization diagram of BA networks due to the positive correlation between node degrees and 
natural frequencies of the oscillations, indicating the evidence of an explosive transition towards 
synchronization of relaxation oscillators system. In contrast to the results in BA networks, in more 
homogeneous ER networks the synchronization transition is always of continuous type regardless of 
the width of the frequency distribution. Moreover, we consider the effect of degree-mixing patterns 
on the nature of the synchronization transition, and find that the degree assortativity is unfavorable 
for the occurrence of such an explosive transition. 

PACS numbers: 89.75.Hc, 05.45.Xt, 89.75.Kd 
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I. INTRODUCTION 

Synchronization is an emerging phenomenon of an en- 
semble of interacting dynamical units that is ubiquitous 
in nature, such as neurons, fireflies or cardiac pacemak- 
ers [H-Q- Inspired by the seminal works of small- world 
networks by Watts and Strogatz^J] and scale-free net- 
works by Barabasi and Albert |5l], synchronization on 
complex networks has been widely studied [6- 8]. These 
studies have revealed that the topology of a network has 
strong influences on the onset of synchronization [9l-[l3|. 
path towards svnchronization and the stability of 
the fully synchronized state However, the con- 

tinuous nature of the synchronization phase transition is 
not affected by the topology of the underlying network, 
even in heterogeneous scale-free networks. 

Recently, explosive transition in complex networks has 
received considerable attention since the discovery of an 
abrupt percolation transition in random networks jl9l.[20| 
and scale- free networks (2ll . [2'2j . although some later 
studies claimed this transition is actually continuous but 
with unusual finite size behavior (23l - [25j . Along this line, 
very recently it was shown that explosive synchroniza- 
tion transitions can take place in scale-free networks of 
phase oscillators [1^ and chaotic oscillators 27]. The 
mechanism responsible for such discontinuous synchro- 
nization transitions is the presence of a positive corre- 
lation between the heterogeneity of the connections and 
the natural frequencies of the oscillators. These studies 



open new perspectives on the research of synchronization 
transitions of other dynamical systems. However, the re- 
search on this topic is still in its infancy and deserves 
more investigations. 

In the present paper, we investigate the criticality of 
phase synchronization transition of coupled FitzHugh- 
Nagumo (FHN) oscillators in Barabasi-Albert (BA) 
scale-free networks Q and Erdos-Renyi (ER) random 
networks [13, [1^. It is shown that the positive corre- 
lation between degrees and natural frequencies of nodes 
can lead to a clear hysteresis loop in synchronization di- 
agram of BA networks and thus signals the occurrence of 
an explosive synchronization transition in heterogeneous 
networks. And a sufficient wide distribution of natural 
frequencies of the oscillators is necessary to induce such 
an explosive transition in BA networks, while in more 
homogenous ER networks the synchronization transition 
is always continuous no matter how wide the distribution 
is. Furthermore, we consider the effect of degree-mixing 
patterns on the nature of the synchronization transition. 
We find that the degree assortativity is unfavorable for 
the occurrence of such an explosive transition. 



II. MODEL 
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Let us consider a network of N coupled non-identical 
FHN oscillators, a representative model of excitable sys- 
tems such as neurons, wherein the dynamics of the ith 
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oscillator is described by the following equations [3C 



N 



ex. 



C'^^Aij{xj—Xi), (1) 
y^ = Xi + ai + ii{t). (2) 



The two dimensionless variables x and y are a voltage- 
like and a recovery-like variable, or in the terminology of 
physical chemistry and semiconductor physics, an activa- 
tor and an inhibitor variable, respectively. The time scale 
ratio £ is much smaller than one (we here set e = 0.01), 
implying x is the fast and y is the slow variable. The 
parameter describes the excitability of the ith unit. 
If \ai\ > 1, the system is excitable, while \ai\ < 1 im- 
plies that the system is oscillatory. The element of ad- 
jacency matrix of the network takes Atj = 1 if nodes i 
and j are connected, and Aij = otherwise. C is the 
coupling constant, and ^i{t) is Gaussian noise that is in- 
dependent for different units and satisfies {£,i{t)) — and 
{^iit)^j{t')) = 2D6^jS{t - t') with noise intensity D. 

To establish a microscopic correlation between the dy- 
namics and the topological properties of nodes, we as- 
sume that the natural frequency uji of the ith unit is an 
increasing function of its degree ki , where ki — Aij 
is the number of nodes that are adjacent to the node. 
Since the oscillation frequency is a decreasing function 
of fli, we thus consider decreases linearly with ki for 
simplicity, i.e.. 



k — k ■ 
tti — 0.99 — 0- ; — 



(3) 



degree in the network, respectively. The factor 6 deter- 
mines the slope of the linear distribution. The larger S is, 
the wider distribution of has, or equivalently, a wider 
frequency distribution. Other parameters N = 200, 
D — 0.005, and the average degree (fc) = 6 are fixed 
in this paper unless otherwise specified. 

To characterize synchronization behavior among the N 
oscillators, we first define the phase of the oscillator i as 

m 



Mt) = 27r- 



27rfc, 



(4) 



where is the time of the fcth firing of the oscillator 
i, which is defined in simulation by threshold crossing of 
Xi{t) = 1.0. Thus, the degree of phase synchronization 



can be measured by calculating r — 

where the vertical bars denote the module and the an- 
gle brackets a temporal averaging. For completely un- 
synchronized motion r ~ 0, while for fully synchronized 
state r ~ 1. 
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FIG. 1; (Color online) Phase synchronization degree r as a 
function of the coupling strength C for different 5 in BA scale- 
free networks (left panels) and ER random networks (right 
panels), (a) and (d) for 5 — 0.1, (b) and (e) for 5 — 0.3, and 
(c) and (f) for S = 0.9. Squares (circles) in (a)-(f) mark the 
forward (backward) simulations, as C is increased (decreased) 
in steps of C = 0.001. Other parameters are A'' = 200, e = 
0.01, D = 0.005, and the average degree (fc) — 6. 



III. RESULTS 

The synchronization diagram is obtained by perform- 
ing both forward and backward simulations. The former 
is done by calculating stationary value of r as varying 
C from to 0.06 in steps of 0.001, and using the final 
configuration of the last simulation run as the initial con- 
dition of the next run, while the latter is performed by 
decreasing C from 0.06 to with the same step. Fig.l 
shows the results of r as a function of C for different val- 
ues of 6 in BA scale-free networks (left panels) and in 
ER random networks (right panels). For BA networks, 
one can find that the nature of synchronization transi- 
tion drastically changes with S. For 6 — 0.1 (Figllja)), 
the results on the forward and backward simulations co- 
incide, implying that the synchronization transition is 
continuous and of a second-order type. Increasing S to 
6 = 0.3 (FiglT](b)), one can see that as C increases r 
abruptly jumps from r^Otor^ilatC^ 0.044, which 
show that a sharp transition takes place for the onset of 
synchronization. On the other hand, the curves corre- 
sponding to the backward simulations also shows a sharp 
transition from the synchronized state to the incoherent 
one at C = 0.021. The two sharp transitions occur at 
different values of C, leading to a strong hysteresis loop 
with respect to the dependence of r on C. Such a fea- 
ture indicates that an explosive first-order synchroniza- 
tion transition arises in BA networks due to the positive 
correlation between degrees and natural frequencies. Fur- 
ther increasing S to S = 0.9 (Figlljc)), the nature of first- 
order phase transition is still present, but the area of the 
hysteresis loop becomes smaller. While for ER random 
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FIG. 2: Area of hysteresis loop in the synchronization di- 
agram of BA scale- free networks as a function of S. Other 
parameters are the same as Fig[T] 



networks (Fig|IJd-f)), the forward and backward simu- 
lations always eoincide regardless of S, and thus such a 
correlation does not induce an explosive synchronization 
transition in ER networks. Therefore, a wide enough dis- 
tribution of natural frequencies and degree heterogeneity 
are both necessary ingredients for the occurrence of such 
an explosive synchronization transition. 

To show how the criticality of synchronization tran- 
sition in BA networks changes with 6, we calculate the 
area of the hysteresis loop in synchronization diagram 
as a function of S, as shown in Figl2l One can see that 
this area equals to zero when 6 < 0.2, implying that the 
synchronization transition is of second-order. When S 
is increased to S = 0.25, the value of this area drasti- 
cally changes to a non-zero value. That is to say, the 
criticality of synchronization transition changes from a 
second-order type to a first-order one at between S = 0.2 
and S — 0.25. With further increasing S, this area show a 
nonmonotonic dependence on S and a maximum area oc- 
curs at (5 = 0.35. This is consistent with the observation 
in FigUl 

To further analyze the change of the order of the syn- 
chronization transition, we have calculated the average 
frequencies ujk of nodes with degree k [k € [fcmin, ^max]) 
along the forward simulation, defined as 



1 



E 

i\ki=} 



2ttT 



(5) 



where Nk is the number of nodes that have degree fc, 
and r ^ 1. By monitoring the variation of Uk with C 
one can clearly observe that how the full synchronization 
state is achieved. In Fig|3Ua) we plot ujk as a function 
of C in BA networks with 6 = 0.3. Before the synchro- 
nization happens, the average frequencies for nodes 
with small degrees and large degrees increase monoton- 
ically, while for those nodes with intermediate degree 
Wfc strongly fluctuate between two transition points. At 
C = 0.044, nodes with any degree class abruptly oscillate 
synchronously, which signals the explosive synchroniza- 
tion shown in FigUfb). For comparison, in Fig|3l^b) and 
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FIG. 3: (Color online) The average frequencies uuk for different 
degrees k along the forward simulation, (a) BA networks: 
S = 0.3; (b) BA networks: 5 = 0.1; (c) ER networks: S = 0.3. 
The arrows in (a-c) indicate the decreasing order of degrees. 
Other parameters are the same as Fig[T] 
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FIG. 4: (Color online) Synchronization diagram in BA scale- 
free networks without any correlation and with the negative 
correlation between degrees and natural frequencies. Other 
parameters are the same as Fig[T] 



FigOl^c) we also plot uJk as a function of C in BA net- 
works with (5 = 0.1 and in ER networks with S = 0.3, 
respectively. One can see that those nodes with large de- 
gree first become locked in frequencies while nodes with 
the small degree classes achieve full synchronization later 
on. This indicates that the synchronization transition is 
continuous. 
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FIG. 5: (Color online) Effect of degree-degree correlations 
on the criticality of the synchronization transition, (a) Syn- 
chronization diagram at S = 0.5 for different degree-mixing 
coefficient rt- (b) Area of hysteresis loop in synchronization 
diagram as a function of 5 for different rk- = 0.0 (black 
squares), = —0.3 (red circles), ~ 0.1 (green triangles), 
and Vk ~ 0.2 (blue diamonds). Other parameters are the 
same as Fig[T] 



Next we will illustrate whether the positive correla- 
tion between degrees and natural frequencies is responsi- 
ble for the explosive synchronization. In Fig|4l we show 
the results of r as a function of C in BA networks, 
where the correlation is destroyed by randomly shuffle 
the values of a^. One can see that the above-mentioned 
explosive synchronization transition disappears and the 
transition becomes a second-order type. Furthermore, 
we consider the case of negative correlation between 
degrees and natural frequencies. To the end, we set 

= 0.69 -I- 0.3{ki - femin)/ (fcmax " ^min) whilc othcr pa- 
rameters keep the same as FiglIJa). The results are also 
shown in FigU] and the transition is continuous. There- 
fore, we can safely conclude that such an explosive syn- 
chronization transition arises due to the positive correla- 
tion between degrees and natural frequencies. 

Lastly, we will consider the effect of degree-degree cor- 
relations on the criticality of the synchronization transi- 
tion. It has been witnessed that many real networks dis- 
play different degree-mixing patterns [32j . To measure 
the degree of the correlation, Newman defined a degree- 
mixing coefficient as [s^l 



rk = 



(6) 



where ji and ki are the degrees of nodes at the two ends 
of the ith link with i — 1 , • • • , M [M is the number of 
total links in the network), = indicates that there is 
no degree correlation, while > (< 0) indicates that 
a network is assortatively (disassortatively) mixed by de- 
gree. A previous study has revealed that degree-mixing 
pattern plays an important role on synchronization [33| . 
To generate different degree-mixing networks, we employ 
an algorithm proposed in 3^. At each elementary step, 
two links in a given network with four different nodes are 
randomly selected. To get an assortative network, the 
links are rewired in such a way that one link connects 
the two nodes with the smaller degrees and the other 
connects the two nodes with the larger degrees. Multiple 



connections are forbidden in this process. Repeat this op- 
eration until an assortative network is generated without 
changing the node degrees of the original network. Simi- 
larly, a disassortative network can be produced with the 
rewiring operation in the mirror method. We start from 
BA scale-free networks with a neutrally degree-mixing 
pattern, and produce some groups of degree-mixing net- 
works by performing the above algorithm. FiglSJa) plots 
the synchronization diagram for different rk at a fixed 
5 = 0.5. One can see that for rk = —0.3 the explosive 
synchronization transition persists. For rk = 0.1 the dis- 
continuous nature of the transition does not change, but 
the area of hysteresis loop become rather small. While 
for Tk = 0.2, the forward and backward simulations coin- 
cide and the synchronization transition becomes contin- 
uous. In Figl5]Jb), we plot the area of hysteresis loop as 
a function of 5 for different rk ■ We find that the case for 
Tfc = —0.3 is similar to that of rk = 0.0. For rk — 0.1 
the area is larger than zero when & > 0.045 but its value 
is very small and has only the order of 10"'^. While for 
Tfc = 0.2, the area is always zero regardless of 6, im- 
plying the nature of the synchronization transition has 
been changed essentially. Therefore, the degree-mixing 
patterns have a significant impact on the criticality of 
the synchronization transition. For a disassortative net- 
work, the explosive nature of the transition appears if S 
is larger than a threshold value, while for a assortative 
network whose degree-mixing coefficient is larger than 
0.2, the explosive transition is absent and the transition 
becomes continuous. Since dissortativity implies nodes 
with larger degrees tend to connect to those nodes with 
smaller degrees, local star configurations are abundant in 
a disassortative network. Thus, it seems to suggest that 
the local star configurations are important for resulting 
in such an explosive phenomenon. 



IV. SUMMARY 

In summary, we have studied the effect of a micro- 
scopic correlation between degrees and natural frequen- 
cies of FHN oscillators on the criticality of synchroniza- 
tion transition in networks of BA and ER models. We 
find that, if the width of frequency distribution is larger 
than a critical value, the positive correlation degrees and 
oscillation frequencies can lead to the first-order syn- 
chronization transition in BA networks. While in more 
homogeneous ER network, such an explosive transition 
dose ont appear no matter how wide the distribution is. 
Therefore, the positive correlation between degrees and 
oscillation frequencies and their heterogeneities are both 
necessary conditions for such an explosive phenomenon. 
Moreover, we have shown the patterns of degree-degree 
correlations have a significant impact on the nature of 
the synchronization transition. In a disassortative net- 
work such an explosive phenomenon persists, while in 
an assortative network the transition becomes continuous 
type if the degree of assortative correlation is relatively 
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larg e. Our results generalizeprevious findings in phase 
[26| and chaotic oscillators [23| to relaxation oscillators 
with time-scale separation, and suggests that the mech- 
anism for generating the discontinuous synchronization 
transition may be universal. In addition, an explosive 
transition may imply that the unsynchronized oscillation 
state can be metastable with respect to the full synchro- 
nized oscillation state near the synchronization transition 
point [35I . The dynamics of spontaneous synchronization 
from the metastable state to the stable full synchronized 
state without any change of system parameters may be 



interesting and deserves further investigations. 
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